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Abstract 

Often, geometries with horizons offer insights into the intricate relation- 
ships between general relativity and quantum physics. However, some sub- 
tle aspects of gravitating quantum systems might be difficult to ascertain 
using static backgrounds, since quantum mechanics incorporates dynamic 
measurability constraints (such as the uncertainty principle, etc.). 

For this reason, the behaviors of quantum systems on a dynamic black 
hole background are explored in this paper. The velocities and trajectories 
of representative outgoing, ingoing, and stationary classical particles are 
calculated and contrasted, and the dynamics of simple quantum fields (both 
massless and massive) on the space-time are examined. Invariant densities 
associated with the quantum fields are exhibited on the Penrose diagram 
that represents the excreting black hole. 

Furthermore, a generic approach for the consistent mutual gravitation 
of quanta in a manner that reproduces the given geometry is developed. 
The dynamics of the mutually gravitating quantum fields are expressed in 
terms of the affine parameter that describes local motions of a given quan- 
tum type on the space-time. Algebraic equations that relate the energy- 
momentum densities of the quantum fields to Einstein's tensor can then 
be developed. An example mutually gravitating system of macroscopically 
coherent quanta along with a core gravitating field is demonstrated. Since 
the approach is generic and algebraic, it can be used to represent a variety 
of systems with specified boundary conditions. 
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1 Introduction 



There have been remarkably few direct experimental observations of the 
quantum behaviors of gravitating systems. Experiments by Overhauser, 
et.al.[lj, have demonstrated that those dynamic gravitational fields local 
to the Earth's surface do not break the coherence of gravitating neutrons, 
which give interference results consistent with the principle of equivalence. 
Those experimental results involve both Newton's gravitational constant 
Gn and Planck's constant h in a single equation form. Such results, along 
with other observed phenomena, require that gravitating sub-clusters can 
maintain quantum coherence while having their internal dynamics influ- 
enced by dis-entangled clusters co-contributing to the local gravity. Al- 
though justifications of the redshift of gravitating photons need not make 
use of Planck's constant, those quanta do maintain their coherence during 
extended interactions with dynamic gravitational fields (e.g. the cosmic 
microwave background, local gravitational redshift measurements, etc.). 
Similarly, interacting gravitating systems continue to gravitate after dis- 
entanglement, which motivates a formulation that incorporates straight- 
forward cluster decomposability within macroscopic gravitational environ- 
ments. 

Generally, the dis-entanglement of relativistic dynamic quantum clus- 
ters that is necessary for correspondence with classical dynamics requires 
that the geometric aspects of the kinematics associated with the Lorentz 
/Poincare transformation properties of a given cluster must be separate 
and distinct from the internal coherent descriptions and off-shell analytic 
behaviors of disparate clusters. The methods utilized in this paper incorpo- 
rate the techniques developed in the establishment of cluster decomposable 
formulations in relativistic few particle scattering [2]. Those formulations 
exhibit the expected dis-entanglement of interacting quantum scattering 
states needed for classical correspondence properties, in spite of the kine- 
matic complications introduced by the non-linear nature of the relativistic 
energy-momentum dispersion relation for massive systems. The solution 
requires proper cluster independence of the geometric parameters describ- 
ing the kinematics between subsystems from the internal quantum dynam- 
ics associated with the description of an interacting system in terms of the 
boundary (i.e., only ^//-interacting) states. Cluster-decomposable rela- 
tivistic scattering theory is most directly realized as follows: 



• the clusters should be characterized using the channel decomposi- 
tion that Faddeev[3] developed for describing cluster-decomposable 
unitary non-relativistic systems; 

• the kinematics between external and intermediate quantum states 
(off- diagonal dynamics) should insure Lorentz frame (3- velocity) con- 
servation rather than 3-momentum conservation. Velocity conserva- 
tion has been referred to as the point form representation of the 
dynamics by Dirac^j and others, whereas momentum conservation 
has been referred to as the contact form representation; 

• the internal dynamics of the intermediate quantum states (or de- 
scriptions of the interacting state in terms of the complete basis of 
boundary states) for the various non-interacting clusters should be 
independent; 

• the complex analytic extension of the dynamic invariant energy of 
a given cluster (the off-shell behavior) should only parametrically 
affect the kinematics of the other clusters, i.e., the internal dynamics 
of one cluster should not alter the energy spectrum of another. 

Formulations of scattering theory with these properties have been shown 
to have appropriate non-relativistic behaviors [5 J as well as to give the ex- 
pected results from perturbative representations such as quantum electro- 
dynamics [5]. 

The relevant quantum dynamics is most conveniently expressed in terms 
of the affine parameter associated with the gravitating particles/fields. The 
affine or proper time derivative of local physical parameters can be ex- 
pressed in terms of the substantive derivative 4^ = u 13 -^ = u^dp. This 
then makes the incorporation of the principles of equivalence/relativity 
straightforward, resulting in substantive quantum flows on the geometry [7j. 

Macroscopic quantum fluids (e.g. superfluids and superconductors) 
maintain persistent quantum flows that satisfy quantization conditions in 
the non-inertial environments of most laboratory measurements. For in- 
stance, vortices of superflow with quantized circulation maintain the an- 
gular momentum of a rotating vessel of liquid helium cooled below the 
superfluid transition temperature. (It would be quite interesting to see if 
the precession of vortices with quantized circulation — in a gravitational 
field Gn is a measurable phenomenon.) Fluid continuity directly follows 



from the equations of motion resulting from Lagrangians constructed using 
substantive derivatives. Superfluid flows then result from the local gauge 
invariance of those Lagrangian forms. Thus, substantive flows have been 
quite useful in describing macroscopic behaviors of quantum fluids. 

Substantive quantum flows inherently articulate the philosophical ap- 
proach of proper time formulations [8 J that presume the primacy of coor- 
dinates associated with the system being influenced in descriptions of the 
interaction. One expects that the affects of external interactions upon a 
given system can be most fundamentally understood in terms of the dynam- 
ical parameters of that impacted system. Thus, when combined with the 
principle of equivalence, the dynamics of free quantum fields in a gravitat- 
ing environment can be most elegantly described using the affine (proper) 
parameters of that quantum field. Furthermore, formulations that exem- 
plify cluster decomposibility and substantive affine flows are most directly 
realized by developing linear spinor fields consistent with those in reference 
[9]. The form of the gravitating quanta developed in Section H] will involve 
affine flows of independent states motivated by the previous discussions. 

The general incorporation of quantum mechanics into gravitating for- 
mulations requires that the descriptions be dynamic, since the relationships 
between temporal and kinematic parameters have inherent constraints due 
to the quantum measurability problem (e.g. the uncetainty principle, non- 
commuting measurables, etc.). Therefore, the background geometry will be 
taken to be that of a dynamic excreting spherically symmetric space-time 
with a metric form given by 



ds 2 = - (l - (dct) 2 +2 ] J^^-dct dr+dr 2 +r 2 {d9 2 +sin 2 6 dcf). 

(1.1) 

In this equation, R M (ct) = 2G^M{ct) / c 2 is a time-dependent form of the 
Schwarzschild radius that will be referred to as the radial mass scale. For 
this dynamic black hole, the radial mass scale is not a light-like surface, 
and is not coincident with the horizon (which is a light-like surface). How- 
ever, a traversing outgoing light-like trajectory is momentarily at rest in 
the radial coordinate r as the radial mass scale moves past that trajectory. 
The trajectory of the horizon is determined by calculating the null radially 
outgoing geodesic ds 2 = that intersects the radial mass scale as the mass 
of the black hole vanishes. This metric was developed as a dynamic exten- 
sion of the river model of (static) black holes discussed in the literature! 10 1. 



The asymptotic (r — > oo) form of the metric is that of a Minkowski space- 
time. Therefore, the coordinates (ct, r, 8, <f>) are the temporal, radial, and 
angular coordinates of an asymptotic observer. 

The temporal parameter t used to describe the dynamics of the black 
hole metric 11.11 is not singular at the horizon. More generally, physical 
curvatures are non-singular [11] away from r = 0. The contravariant form 
for the Einstein tensor calculated from this metric can be expressed in terms 
of the excretion rate R,M(ct) and the dimensionless parameter ( = Rm ^ ; 
taking the form 
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Since several components of the Einstein tensor are non-vanishing, this 
dynamic geometry provides an excellent laboratory for the exploration of 
self-consistent gravitating phenomena. 

The subsequent developments will exhibit the behaviors of classical and 
quantum objects on the black hole background described in Eqns. II .11 and 
11.21 Section [2] will exhibit solutions to the geodesic equations for the ge- 
ometry, as well as demonstrate typical trajectories of massive outgoing, 
ingoing, and stationary classical particles on the Penrose diagram of the 
space-time. The expected behaviors of classical massless particles have 
been confirmed in a previous paper [13]. The behavior of a spherically 
symmetric Klein-Gordon field will be exhibited in Section [3j The form of 
the Klein-Gordon field will be found to be clumsy for constructing mutu- 
ally gravitating fields, so a more convenient set of cluster-decomposable 
scalar fields will be developed in Section |H Finally, in Section an ex- 
ample mutually gravitating macroscopic system of quanta is exhibited and 
discussed. 



2 Geodesic Motion 



2.1 Four velocities 

Geodesic motion is at the foundation of the principle of equivalence. 
The geodesic equation describes the evolution of four-velocities on the ge- 
ometry. In all calculations, four- velocities of both massive (^f—) and mass- 
less (^fr-) particles will be taken to be dimensionless. The form satisfied 
by the radial component of four- velocities on the metric 11.11 is given by 



r = -C /2 u'± J(u°y-e m (2.1) 



u 

where Q m = j ^ ^ . In particular, stationary massive particles sat- 
isfy u° = l,u r = — C 1 / 2 - Those four-velocities associated with geodesic 
motion satisfy the geodesic equations: 

dv° C 1 / 2 



d\ 2r 

and 



(C 1/ V + « r ) a > (2-2) 



Away from the horizon ( = (h oi the ingoing horizon ( = (i (which 
of course define the geodesic of a particular outgoing or ingoing massless 
particle, respectively), the four- velocities of freely gravitating particles can 
be reparameterized in terms of the dimensionless variable ( in the black 
hole geometry. This reparameterization involves the transformation 

d( = Rmu° - (u r ^ ^ 

dX r 

Generally, affine derivatives can be expressed in terms of the substantive 
(or flow) derivative 4- = u^dp, which, using Eqns. 12.21 and 12.31 implies 
that gravitating four-velocities can be expressed as functions of the single 
parameter (, i.e. u@ = u^(Q. 

Therefore, numerical solutions of the geodesic equations can be di- 
rectly obtained as ordinary differential equations. For all numerical cal- 
culations, the (dimensionless) excretion rate will be taken to be constant 



Rm = 0.1. This produces an outgoing horizon dressing the singularity 
r = at Ch — 1-177 and an ingoing horizon of outermost communication 
with the singularity at Q = 0.078. Plots of the four-velocities of massless 
particles with asymptotic values u°(0) = 1 are demonstrated in Figured! 
while plots of the four-velocities of massive particles with asymptotic values 
u°(0) = 1.2 are demonstrated in Figure [2j The four- velocities of outgoing 
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Figure 1: Four velocities for massless particles as a function of (. Upper 
curves represent u° and lower u T 




Figure 2: Four velocities for massive particles as a function of (. Upper 
curves represent u° and lower u r 

particles are seen to have singular behavior at the horizon r = i?#(c£), 
while the ingoing particles have singular behavior at the ingoing horizon 
r = Rj(ct). The radial components of the four- velocity of outgoing parti- 
cles are seen to change sign at the radial mass scale (m = 1- Stationary 
massive particles are described by a vanishing value for the radical in Eqn. 
12.31 The plot of stationary particles in the space-time is demonstrated 
in Figure [3j Any of the above-calculated forms for the four-velocities of 




Figure 3: Stationary massive u° (upper plot) and u r (lower plot). 

massless and massive systems can then be used to develop the classical 
trajectories of gravitating particles. 

2.2 Form of conformal coordinates 

Space-time diagrams are always useful in the visualization of dynamics 
in a given geometry. In particular, Penrose diagrams are convenient for 
examining the large-scale causal structure of a geometry because of the 
following properties: 

... Penrose diagrams map the entire space-time onto a single finite page, 
and 

... Penrose diagrams preserve the slope of light-like trajectories at ±unity. 

Since light-like curves are easy to plot on a Penrose diagram, the causal 
structure of the geometry can be directly observed from the diagram, and 
potential causal relationships between locations can be immediately ascer- 
tained. In order to construct the Penrose diagram for a given space-time, 
conformal coordinates, which give unity slope for light-like curves, must 
be found. For Minkowski space-time, the coordinates (ct, r) are already 
conformal. However, for the metric 11.11 conformal coordinates must be 
constructed, since clearly null geodesies will not have unity slope for those 
coordinates. 



For a black hole satisfying Eqn. 11.11 with a constant rate of mass accre- 
tion/excretion Rm = 0, the conformal coordinates can be determined in a 
straightforward manner. The forms of the conformal temporal and radial 
coordinates that are used for the construction of the Penrose diagram have 
been developed in a companion paper [12], and are given by 
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(2.5) 

This equation relates the space-time coordinates of an asymptotic observer 
(ct,r) with the conformal coordinates (c£*,r*). The Penrose diagram of a 
spatially coherent black hole that excretes at a constant rate Rm until it 
vanishes at t = 0, is shown in Figure HI The singularity of the black hole 
is represented by the bounding space-like curve r = on the left-hand 
side of the diagram. At t = 0, this singularity vanishes, and the curve 
r = becomes time-like in the final Minkowski space-time (which is rep- 
resented by the upper right quadrant of the diagram). No communication 
to the left of the light-like horizon Rh can escape hitting the singular- 
ity. Likewise, no communication to the right of the light-like incoming 
horizon Rj can communicate with the singularity. The radial mass scale 
RM^ct) = 2G]<qM{ct) I c 2 is seen to itself be a time-like trajectory where out- 
going light-like trajectories remain temporarily stationary in coordinate r. 
The conformal coordinates (ct*, r*) in Eqn. 12.51 are chosen to correspond to 
the Minkowski coordinates (ct, r) (which are, of course, also conformal) of 
the asymptotic observers far from the black hole. Therefore, the bounding 
light-like curves on the right-hand quadrants of the diagram correspond to 
the Minkowski space-time of distant observers. 



2.3 Classical trajectories 

The classical trajectories of freely falling bodies can be directly repre- 
sented on the Penrose diagram by integrating the four-velocities subject 
to the appropriate initial conditions. All light-like trajectories are indeed 
found to always have slopes of ±unity[13j. Massive particle trajectories 



Figure 4: Penrose diagram for a black hole that evaporates steadily to zero 
mass at ct — 0. Red curves (running vertically in the right hand region) 
represent curves of constant r successively in tenths of units, units of length, 
and decades of units. The green curves represent curves of constant ct in 
units of length. The horizontal solid blue line on the right represents the 
end of excretion ct = (Co = 0), while the horizontal red curve on the 
left represents the space-like singularity of the black hole r = (( s = oo). 
The dashed blue line labeled Rh ((h ~ 1.177) represents the horizon, that 
labeled Rm (Cm — 1) represents the radial mass scale, and that labeled Ri 
(C/ ~ 0.078) represents the ingoing horizon. 



are demonstrated in Figure [5j In the first diagram, various outgoing and 




Figure 5: Classical trajectories for massive gravitating particles. The 
Penrose diagram on the left demonstrates outgoing trajectories (a) and 

(b) terminating on the singularity at ct = —5 , outgoing trajectories (c) 
and (d) passing through r = +5 at ct = 0, and ingoing trajectories (e) 
terminating at ct = —5, and (f) passing through r = +5 at ct = 0. The 
diagram on the right demonstrates stationary trajectories terminating on 
the singularity at ct = —5, and passing through r = +5 at ct = 0. 

ingoing trajectories are plotted. Trajectory (a) is an outgoing mass that ul- 
timately hits the singularity, while trajectory (b) is an outgoing mass that 
hits the singularity at the same time, but moving more slowly. Trajectory 

(c) is an outgoing trajectory that escapes hitting the singularity, while tra- 
jectory (d) is a faster outgoing trajectory that passes through the point 
(ct = 0, r = 5) at the same time. Trajectory (e) is an ingoing mass that 
hits the singularity at the same time as trajectories (a) and (b). Trajectory 
(f) is an ingoing mass that is external to the ingoing horizon Ri, and passes 
through the point (ct = 0, r = 5) at the same time as trajectories (c) and 

(d) . No numerically stable solution for massive outgoing trajectories in 
the region between the horizon and the radial mass scale was found. The 
second diagram plots the trajectories of stationary masses. The trajectory 
on the left of the second diagram represents a stationary mass that hits the 
singularity at the same time as trajectories (a), (b), and (e) on the first 



diagram. The trajectory on the right of the second diagram represents 
a stationary mass that passes through the point (ct = 0,r = 5) at the 
same time as trajectories (c), (d), and (f) on the first diagram. All tra- 
jectories that escape the singularity are seen to smoothly transition from 
gravitating projectiles to inertial free motion as the space-time transitions 
from dynamic black hole to Minkowski space-time at t = 0. It is interest- 
ing to note that a stationary particle which remains in the exterior region 
of the black hole geometry apparently follows the same trajectory using 
conformal coordinates as it would in Minkowski space-time. 



3 Klein-Gordon Field 



The primary focus of this paper is the examination of quantum behav- 
iors on the dynamic space-time defined by the metric 11.11 It is instructive 
to first examine the quantum mechanics of a well-understood scalar field 
on the space-time background described by this metric. The Lagrangian 
of a non-interacting gravitating massless scalar Klein-Gordon field x will 
take the form 

1 



^-g g^d^x f {x)d vX {x). 



(3-1) 



The massless field is chosen so that no additional scale is introduced into 
the problem by including a mass for the quanta. As shown in reference 
j , one can perform a partial wave expansion on \ 



(3.2) 



to obtain equations [15] describing the dynamics of the field ipi 
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(3-3) 

S-wave (JL = 0) solutions of Eqn. 13.31 have been numerically calculated 
for the excreting black hole. Boundary conditions at ct = require that the 
field i/jq and its derivatives must appropriately match those of a free-space 



Minkowski Klein-Gordon field near that region. The matching Minkowski 
space solution for the field away from r = ct is given by 



ip M (ct,r) 



A + Blog 



r + ct 
r — ct 



Q(r - ct). 



(3.4) 



Plots of s-wave solutions to Eqn. 13.31 on the black hole background are 
shown in Figure The field was chosen to have a very small (though 
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Figure 6: Plots of ^(C) and ^'(C) f° r a uiassless Klein-Gordon field on the 
dynamic black hole background, where ( = RM{ct)/r. Field normalization 
units are arbitrary. 

non- vanishing) value asymptotically (( = 0). The wavefunction is contin- 
uous throughout the space-time, but its derivative is seen to have singular 
behavior near the ingoing horizon Rj (£j pa 0.078) and the horizon Ru 
((h ~ 1.177). The magnitude of the scalar field is seen to be largest within 
and near the horizon. 

To gain physical insight into the distribution of the scalar field on 
the global space-time, a density plot on the Penrose diagram is instruc- 
tive. For generating such a plot, the computer was instructed to place a 
pixel at a point on the Penrose diagram corresponding to the asymptotic 
observer's coordinate label (ct, r) if the density multiplied by a random 
number bounded by and 1 was larger than the value calculated at a 
normalization point. The plot then demonstrates relative measures of the 
density throughout the space-time. The quantity plotted in Figure [7] is 
the probability density \x(ct,r)\ 2 . In the first density plot, the normal- 
ization point (represented by the center of the small circle) was placed on 
the radial mass scale curve Ruict = —5). The second curve represents the 
same density, only placing the normalization point just outside of the inner 
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Figure 7: Penrose density plot of a spherically symmetric gravitating 
massless Klein Gordon field (left) and that same field's correspondence 
with the Minkowski space-time solution (right). 

horizon, which allows smaller values of the field density to be plotted. The 
weak field black hole solution for the field is seen to reflect the Minkowski 
space-time solution that matches that field on the line (volume) ct — 0. 

An intriguing question is whether the spherically symmetric massless 
Klein-Gordon field can self-gravitate to produce this black hole geome- 
try. The energy-momentum tensor generated by the Lagrangian 13.31 can 
be calculated several ways. For this field, the action was calculated from 
the Lagrangian using W = J Cd 4 x . The energy-momentum tensor was 
then calculated by examining the behavior of this action under variations 
of the metric, 5W = hj d 4 x^—g T^ u 5g^ u . One can then compare this 
tensor to the Einstein tensor 11.21 to examine how the scalar field might 
contribute to the gravitating system. In particular, it is convenient to ex- 
amine the physical scalar density invariant g^ v T^ u (i.e. the trace of the 
scalar field's energy- momentum tensor), which would be directly related 
to the Ricci scalar 1Z of the space-time if the Klein-Gordon field was self- 
gravitating. The Klein-Gordon density invariant and the Ricci scalar are 
plotted in Figure [SJ The density invariant is singular at the horizon and 
ingoing horizon, while the Ricci scalar is not. The Klein-Gordon field 
here calculated was not expected to directly construct a self-gravitating 
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Figure 8: Plots of invariant scalar densities associated with gravitating 
quanta on the black hole background. The left diagram plots r 4 times 
the trace of the energy-momentum tensor of the Klein-Gordon field. The 
diagram on the right plots r 2 times the Ricci scalar of the dynamic black 
hole. The horizontal axis for both diagrams is the dimensionless parameter 
C = R M (ct)/r. 



quantum consistent solution to Einstein's equation. The non-linear form 
of the Klein-Gordon energy-momentum tensor in derivatives makes super- 
positions clumsy. That form also makes the dependency of the density 
invariant upon the radial coordinate r irreparably different from that of 
the Ricci scalar. In addition, the form does not intuitively incorporate 
some of the properties that one would expect of a self-gravitating field. In 
the next section, a foundation will be developed for the exploration of grav- 
itating quantum fields on a space-time background that can be superposed 
in a manner that might construct mutually gravitating objects. 

4 Cluster-Decomposable Gravitating Scalar 
Fields 

Useful gravitating quanta that mutually contribute to a particular space- 
time configuration are expected to satisfy certain physical properties: 

... it should be possible to have dis-entangled components that gravitate, 
and can co-contribute to a mutually gravitating environment; 

... the energy densities of the components should linearly contribute to 
the overall energy density used to drive the Einstein equation for the 



gravitating environment; 



... there should be components that can maintain quantum coherence, 
satisfying probability flux conservation and expected energy/time- 
momentum/position phase relationships; 

... gravitating quantum systems should be able to maintain observed 
linearity as previously dis-entangled components interact, forming 
changing configurations of those systems; 

... the collective gravitational field should incorporate the principle of 
equivalence with regards to the gravitation of a given component, 
i.e., one should obtain expected behaviors when performing scattering 
and spectroscopic experiments in a locally inertial freely falling frame. 

Interacting few-particle quantum systems exhibit peculiar behaviors, 
such as the Efimov effect [16] (which demonstrates long range coherence 
for short ranged interactions), even for weakly interacting non-relativistic 
systems. The non-linear behavior of the energy-momentum relationship in 
Minkowski space-time further complicates attempts to explicitly demon- 
strate how interacting quantum systems can dis-entangle. The problem 
is that the complicated relativistic kinematics between entangled clusters 
will generally alter descriptions of the dynamic properties of one clus- 
ter when another cluster is interacting. For instance, in naive relativis- 
tic formulations an electron interacting on the moon would substantially 
affect the energy spectrum of a hydrogen atom here on earth. How- 
ever, cluster-decomposable formulations of quantum scattering maintain 
relativistic covariance, and incorporate correspondence with the classical 
dis-entanglement of clusters^ [5]. The flat space-time solution of cluster- 
decomposable relativistic quanta requires that one properly separate the 
geometry (kinematics) from the off-shell quantum dynamics. As was men- 
tioned in the Introduction, this is done by formulating the dynamics while 
incorporating the following properties: 

(i) a recognition of the difference between off-shell and off-diagonal 
descriptions of intermediate states. The off-shell behavior describes 
the analytic structure of the invariant energy description of the scat- 
tering amplitudes in the complex plane. Bound states manifest as 
poles in the off-shell parameter. The off-diagonal behavior describes 



how the fully interacting system can be described in terms of a com- 
plete set of boundary states. The intermediate (virtual) state dy- 
namics is separate from the off-energy-shell parametric behavior; 

(ii) the off-diagonal description of intermediate quantum states should 
use Lorentz frame conservation (parameterized by preserving the 
three components of the four-velocity u, or so called point form dy- 
namics) instead of momentum p conservation (or contact form dy- 
namics); 

(iii) the formulation should use parametric (geometric rather than 
quantum dynamic) descriptors of cluster kinematics, in terms of the 
invariant energies u ■ P for the overall system, and u a ■ P a for cluster 
a. 

In order to most directly infuse these properties into the formulation, quan- 
tum states with explicit linear behaviors in derivatives will be constructed. 



4.1 Dynamical equations of gravitating quanta 

As previously mentioned, Lagrangians that utilize substantive derivatives [7] 
most directly incorporate the spirit of the principle of equivalence. The 
chosen form for the Lagrangian of a gravitating non-interacting cluster a 
in the geometry will be taken to be 

£ a = L a = -y/^g {—uP^ld^a - (dpipDipa) - m a c 2 ip* a ijj a ]. (4.1) 

Such a Lagrangian form can be directly generalized to incorporate linear 
spinor fields [9]. 

Writing the functions ip a as complex parameters ipa = \^a\ e i t ne 
Euler-Lagrange equations then insure probability conservation 

1 =5 /J (V=ff|V'«|V) = ) (4.2) 



and phase coherence 

u%i a = -^f. (4.3) 

Therefore, systems described by the Lagrangian form Eqn. 14.11 satisfy 
the dynamics expected of inertial quantum fields. Substituting the phase 



coherence equation 14.31 back into the original form for the Lagrangian, it 
can be shown that the extremum Lagrangian for the gravitating quanta 
vanishes C a = 0. 

If one defines E m as the proper value (standard form invariant of the 
little group [TT] ) of the energy of the particle (i.e. mc 2 for massive particles, 
E for massless particles), a general solution to Eqn. 14.31 is given by 



he he 



'l-()u° + (J( + Q)u 



Ort = ^(u r ~Qu°) = ^ Uy/C ~ Q)u° + U ' 

he he L 

where the covariant forms of the 4- velocities are given by 

Uo = (C - l)u° + C 1/2 ' 



(4.4) 
(4.5) 



(4.6) 



The four- velocities satisfy the usual normalization u^up = —Q m . 

The integrability of a function of several variables means that the func- 
tion is in fact well described in terms of those variables, satisfying the 
analytic property of equality of second partial derivatives regardless of or- 
der. The integrability of the phase [d r ,d ]^(ct,r) = implies that the 
function Q must satisfy 



^dpQ = d u r - d r u - (d u ) Q. 



(4.7) 



Since the four-velocities depend only on the the dimensionless parameter 



C, uf 



u 



1 Rm (ct) ■ 



the functions Q = Q(C) can likewise be shown to 



depend only upon the dimensionless parameter (. For the metric form II. 1\ 
derivatives of the geodesic four- velocities (from Eqns. 12.21 and 12.31) satisfy 



d a u 



and 



e m c 1/2 ^m + C 3/2 

r 2 R M u°-(u r 



d u r — d r u = 0. 



(4.8) 



(4.9) 



Numerical solutions for the magnitudes and phases of the fields can be 
developed by defining the dimensionless parameters ^(C)! 2 by the equation 



\ip(ct, r)| 2 = , and the reduced phases (which carry the dimension 
of inverse Compton wavelength) by the equation £(ct,r) = r£(£). It is 
interesting to note that ^(Ol 2 an< ^ Q(0 have the same form for their 
differential equations. 

For massless particles, one can immediately determine that the parame- 
ters |^>(C)| 2 and Q(C) are constants. Analytic solutions can also be obtained 
for stationary particles. For stationary massive particles, the solutions take 
the form 

l^(C)l 2 = ( Rm + C 3 J 2 ) 1/3 = QsiO u 10) 

\MCo)\ 2 \R M + C 3/2 J QsiO' 

Numerical forms describing the quantum dynamics of outgoing and ingoing 
massive quanta can directly be obtained by solving the straightforward 
ordinary differential equations for the phase factor Q(C) (as well as the 
reduced squared wavefunction ^(Ol 2 ) using Eqn. 14.71 and integrating 
Eqn. 14.41 to obtain the phase £. 



4.2 Energy- momentum tensors of gravitating quanta 

The form of the contribution of the gravitating quanta to the energy- 
momentum tensor can be directly calculated using the standard Legendre 
transformation from the Lagrangian form to the Hamiltonian form of the 
dynamics of a multi-component field \'- 

T^ = Z^X-%C. (4.1!) 

The Euler-Lagrange equations imply that the extremal form of the La- 
grangian 14.11 vanishes C a = for any quantum a. Therefore, the energy- 
momentum tensor of quantum a satisfies 

= V^T a \ = V^(hc) u P(d^ a )\ip a \ 2 . (4.12) 

The energy-momentum contribution of quantum a to the overall system is 
seen to depend linearly upon derivatives of the phase, directly linking the 
local phase of the wavefunction a to its contribution to the local energy- 
momentum of the system. 

The form of these generic quanta need not locally conserve energy- 
momentum due to the background gravitation. The generic form for the 



divergence of the of the energy-momentum tensor of a field ip a is given by 



T ati-,(3 = a/ a 1 + complex conjugate. (4.13) 



Neither is the form of this tensor necessarily symmetric, satisfying 

T a 0r -T a r0 = e ma ^\^\ 2 Q a . (4.14) 

Since the Einstein tensor 11.21 is both geometrically conserved and symmet- 
ric, the collective form of the energy-momentum tensor driving the Einstein 
equation must likewise be locally conserved and symmetric. 



5 Macroscopic Self-Gravitating Quanta 

Although the form of microscopic quantum gravity remains uncertain, 
the quantum mechanics of systems co-gravitating with macroscopic me- 
dia has been phenomenologically explored by common experience as well 
as experiment pj. As a complement to attempts to self-consistently de- 
scribe the microscopic behaviors between mutually gravitating quanta, one 
should be able to gain insight into the fundamentals of gravitation by de- 
veloping micro-physical behaviors that consistently reproduce an example 
dynamic space-time. Therefore the form of mutually gravitating coherent 
sub-clusters that are consistent with the metric Eqn. 11.11 will be developed 
in this section. 

The geometrodynamics of the space-time must consistently co-mingle 
with the micro-physical behaviors for the construction to be viable. The 
Ricci scalar 1Z used in the calculation of the Einstein tensor is directly 
related to a physical invariant given by the trace of the energy-momentum 
tensor T^p. For the space-time metric [TTTl it is given by 



which is non-singular away from a physical singularity at the origin, and 
vanishes for a static radial mass scale Rm —* for finite radial mass scales 
Rm- However, one should note that although the metric Eqn. 11.11 takes 



the same form as a Minkowski space-time as the radial mass scale van- 
ishes, some physical attributes of a black hole that excretes at a steady 
rate become singular when Rm — > in Eqn. 15.11 This paper is not con- 
cerned with how the excretion rate Rm turns off as the radial mass scale 
vanishes, which is what must occur in an actual physical process. Reason- 
able calculations show that one can obtain viable quantum solutions with 
stepwise constant rate equations for radial mass scales down to the Planck 
length. Therefore, gravitating quantum solutions will be examined for the 
constantly excreting space-time represented in the lower quadrants of the 
Penrose diagram in Figure H] for times t < 0. 



5.1 Core quantum field 

The Lagrangian densities for the cluster-decomposable gravitating scalars 
developed in Section H] vanishes for extremal fields. The space-time de- 
scribed by the metric [1J] is spherically symmetric, so that all velocities have 
vanishing angular components u e = = vr. This means that quanta with 
energy-momentum tensors of the form Eqn. 14.121 cannot contribute to the 
non-vanishing components of the Einstein tensor G t \ and G^^. Therefore, 
the vanishing of the Lagrangian densities C a = for the non-interacting 
gravitating quanta, requires the addition of an interacting form for a core 
gravitating field in order to be consistent with Einstein's equation. 

The core gravitating field is expected to relate closely to classical geo- 
metric parameters, and will be represented by real fields with a Lagrangian 
of the form 



(5.2) 



The core fields ip c ± satisfy the same Euler Lagrange equations, but will 
have opposing signs for their energy densities. The core gravitating fields 
are seen to have an interaction term dependent on the local form of the 
metric. 

For the real fields ip c ±, the Euler-Lagrange equations require that 

^^±8^ = 0. (5.3) 



From Eqn. I4.8[ one observes that any massless core fields will satisfy these 
Euler-Lagrange equations. Light-like propagation properties for the core 



gravitating field are, of course, intuitively gratifying. This gives consider- 
able flexibility in choosing forms for the core fields ip c ± that are consistent 
with the requirements of Einstein's equation for this geometry. 
The form of the core field energy-momentum tensor is given by 



T p = -tic 

-'-core fi 



ip c+ u /3 d II il) c+ - ipc-uPd^ipc-} - S^CcorJy/^g. (5.4) 



Einstein's equation then requires that the physical content of the system 
satisfies 

G\ = -87r^|(T core ^ + T ro / M ) (5.5) 



where the Planck length is defined by Lp = -^— c = y^§F and the tensor 
of radiating quanta T ra d is the collective form of the cluster-decomposable 
quanta developed in Section HI As previously stated, since the radiations 
defined in Eqn. 14.11 cannot contribute to G 2 2 = G 3 % for a spherically 
symmetric system (since u 9 = = in the equation 14. 121 for the radiating 
components' contributions), Einstein's equation constrains the form of the 
core gravitating field from Eqn. 11.21 

7? ' W - G 3 3 = -87T^T core 3 3 = Stt^/W^, (5.6) 



4r 2 v/C a ""he™* he 
which yields the equation satisfied by the core gravitating fields 



A solution of Eqn. 15.71 with positive semi-definite probability densities 
can be found. The solution required (i) outgoing massless core quanta, 
(ii) ip c+ = for r < Rn(ct) and (iii) ipe- = for r > Rn{ct). If de- 
scribed in terms of a single core field, the square of that single field simply 
changes signs across the horizon. Solutions are demonstrated for reduced 
core probability density |Y" C (C)| 2 (where \ip c (ct,r)\ 2 = 1 ) and reduced 

core energy density j- T e ° ° e in Figure M The solution for the square 
of the core wavefunction is everywhere non-negative, but vanishes at the 
horizon Ch — 1.177. The core field energy density is negative in regions 
distant from the black hole (for £ < 0.2), positive in the region outside but 
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Figure 9: Reduced wavefunction squared and energy density of core field 
as a function of £ — i?A/(ct)/r. 

near the horizon and again negative inside the horizon. Since the com- 
ponent G 00 vanishes, the sign taken by the energy density of the radiating 
quanta from Section H] must cancel the energy density contributed by the 
core field. It was therefore crucial that the asymptotic form of the core 
field have this negative energy density so that normal gravitating particles 
far from the black hole have the expected sign for their energy densities. 

As was seen for the Klein-Gordon field in Section [3] it is illustrative to 
examine a density plot of the physical densities calculated in Figure [H]on the 
global space-time represented by the Penrose diagram. As was previously 
discussed, only that portion of the Penrose diagram satisfying t < will 
be physically populated. The square of the core field wavefunction and 
energy density are plotted in Figure [TU1 As with all density plots, these plot 
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Figure 10: Penrose density plots of core field probability density and 
energy density. 

represent the relative magnitude of the given density in terms of the relative 
number of pixel points plotted in a given vicinity. The normalization point 



is represent by the small circle in the middle of the diagram near the radial 
mass scale Cm = 1- Since the physical parameters being plotted might take 
on negative values, a blue pixel is plotted if the density is positive relative 
to the absolute value of the normalization scale, while a red pixel is plotted 
if it is negative relative to the negative absolute value of the normalization 
scale. The core field energy density is seen to take on negative values in 
the interior region, as well as small negative values in far regions on the 
right of the Penrose diagram. 

Solutions of Eqn. 15.71 involving ingoing massless core quanta were like- 
wise calculated. The solution for ingoing massless core quanta obtained 
by the author had negative probability densities in a small region near 
the radial mass scale (within numerical accuracy), and were not further 
explored. 

5.2 A mutually-gravitating macroscopic system 

Once the core gravitating field of the excreting black hole has been 
determined from Eqn. I5.7[ the linear behavior of the radiation fields can 
be exploited to construct a mutually/self gravitating system. Substituting 
solutions for the core field calculated from Eqn. 15.71 into the core energy- 
momentum in Eqn. 15.44 this energy-momentum tensor then placed into 
Einstein's equation 15.51 defines a form that must be satisfied by the collec- 
tion of radiating gravitating quanta. 

From Eqn. 14.121 a generic form for the macroscopic gravitating quanta 
is given by 

T rad \ = Y,^c)u^ a )m 2 . (5.8) 

a 

It is very gratifying that in the exterior region of the black hole, there are 
many, many solutions satisfying Einstein's equation [53] due to the algebraic 
and dis-entangled nature of the contributions from the radiating quanta. 
This means that the phenomenologically observed behaviors of dynamic 
spherically symmetric stars, planets, etc. can likely be incorporated using 
mutually gravitating particles as described by the linearly additive Eqn. 

EE 

An example solution for a mutually gravitating macroscopically coher- 
ent system of quanta will be directly demonstrated. Since there are 4 com- 
ponents of the Einstein tensor Eqn. I1.2l that remain to be satisfied by Eqn. 



I5.5[ a minimal algebraic solution should be obtained using 4 particle types 
in Eqn. 15.81 The geometry should have coherent emissions/absorptions of 
otherwise conserved radiating quanta locally contributing to the energy- 
momentum of the system in the form 



T ra / M = (He) [u?(0„£O Wil 2 + 4(d,i 2 )N 2 \^ 2 \ 2 + 



(5.9) 



where N a represents the number of quanta of type a present in the grav- 
itating system. As long as the fields are linearly independent, solutions 
to Einstein's equation 15.51 using macroscopic quantum fields of the form 
in Eqn. 15.91 can always be found. However, one is not guaranteed that 
physical boundary conditions with positive semi-definite number densities 
can be found. It was found to be quite straightforward to obtain physi- 
cally meaningful solutions with positive number densities and positive mass 
quanta in the region exterior to the radial mass scale r > Rm- However, the 
author has yet to develop solutions that have both positive semi-definite 
number densities and positive semi-definite masses in all regions. This is 
likely due to the behavior of the component u r as seen in Eqn. 12.31 since it 
can change its sign in the interior region of the black hole, and therefore af- 
fect the algebraic solutions. An example macroscopic mutually gravitating 
excreting black hole is given below. 

The solution demonstrated in Figures [TT], [TJl dSJ and O involve a 
system with four particle types. Particle types 1 and 2 are stationary 
gravitating quanta with mass m and linearly independent phases demon- 
strated in Fig. [TTJ The reduced phases = — ^' r ^ of the station- 
ary particles are represented in the first two diagrams in the figure. The 
magnitude squared of the reduced wavefunction of the stationary parti- 
cles |^> a (C)| 2 = r2 Lp\ip a (ct, r)| 2 is represented in the third diagram in the 
figure. The singular behavior in the stationary form of |V>(C)| 2 occurs at 
Css — (—Rm) 2 ^ 3 — 0.215 for the excretion rate chosen for this paper, 
which is in the region of the space-time between the radial mass scale Rm 
(Cm = 1) and the ingoing horizon Ri (£j = 0.078). In Figure fT2l particle 
type 3 is an outgoing massless quantum represented by the first diagram 
and particle type 4 is an ingoing massless quantum represented by the sec- 
ond diagram. The reduced wavefunction for massless quanta is constant, 
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Figure 11: Phases of stationary mass 1 and mass 2, and the square of 
reduced stationary wavefunction. 
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Figure 12: Outgoing and ingoing massless particle phases. 



and the reduced phase is expressed in terms of the quantum's asymptotic 
energy E using |(C) = jr^f 2 ^. 

Algebraic solutions for the reduced number densities of the four par- 
ticle types M a {C) = ^rLp r 2 N a (ct, r)\ip a (ct, r)\ 2 (where again E m is the 
proper or affine energy form of a quantum of mass m), are demonstrated 
in Figure O The plot representing a particular particle type is directly 
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Figure 13: Reduced number densities of mutually gravitating quanta. Par- 
ticle types (1) and (2) are stationary massive quanta. Particle type (3) is an 
outgoing massless quantum, while particle type (4) is an ingoing massless 
quantum. The horizontal axis is the dimensionless scale ( = RM{ct)/r 

labeled on the diagram. The affine energy form E ma (associated with the 
affine parameter describing the particle motion) appears in the phase of the 
wavefunction of particle type a as is the case with standard quantum me- 
chanics. For massive states it represents the "co-moving" energy mc 2 , while 
for massless states is represents the "standard state" energy-momentum E 
of the particle observed in the boundary inertial reference frame. One can 
see that particle types 2, 3, and 4 must take on negative values of E m in 
the region just inside of the radial mass scale r = R,M{ct) ((m = 1), and 



outside of the horizon r = Rji{cb) (Ch ~ 1-2), if the quantum number den- 
sities of those particle types are to remain positive semi-definite. Thus, one 
expects to find some negative mass stationary quanta of type 2 just inside 
the radial mass scale, but outside of the horizon. Particle type 1 is seen 
to require a negative mass value inside the horizon near the singularity. 
Since the solutions to Eqn. 15.51 are algebraic and only depend upon Ein- 
stein's equation, the number densities themselves are independent of the 
normalization of the wavefunctions. Only the actual numbers of quanta 
satisfying a particular flux normalization depend upon the form of that 
normalization. Since the problem examined here has unbounded past tem- 
poral extent, normalization properties will be examined using models for a 
complete finite duration life-cycle (accretion, then evaporation) of a black 
hole in subsequent papers. 

It is particularly illustrative to examine density plots of the various 
number densities on the Penrose diagram. Density plots of the square of 
the wavefunction and number density of particle type 1 are displayed in 
Figure [TH The density plots show densities relative to a normalization 
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Figure 14: Penrose density plots of probability density and number density 
x mass of stationary mass 1 

point in the center of the diagram just outside of the horizon indicated by 
the small circle. Once again, blue pixels indicate positive values of that 
density, while red pixels indicate negative values. All probability densities 
are found to be positive semi-definite throughout the global space-time. 
One might note that as previously stated, a positive semi-definite number 
density requires negative mass solutions near the singularity, indicated by 
the red values on the second density plot. 

Scenarios involving various combinations of ingoing, outgoing, and sta- 
tionary massive quanta co-gravitating with ingoing and outgoing massless 



quanta have been calculated. The form of the four-velocities for the given 
geometry are fixed by the particle type (e.g. massless vs. massive). The 
algebraic nature of the quantum solutions make the numerical calculations 
straightforward. To obtain a given solution, one needs to adjust boundary 
values for the phases and velocities in order to match the given physical 
boundary conditions. It is expected that one should be able to describe 
physically mixed solutions (e.g. thermal systems or incoherent sums of so- 
lutions consistent with the boundary conditions) in a direct manner using 
this formulation. 

6 Conclusions 

The behaviors of classical particles and simple quantum fields on a 
dynamic black hole have been both analytically and numerically exam- 
ined. The metric form taken to describe the dynamic black hole directly 
corresponds with a Minkowski space-time asymptotically, and utilizes the 
asymptotic observer's time coordinate to describe the dynamics of the black 
hole without introducing physically singular behavior at the horizon. This 
metric then defines the form for outgoing, ingoing, and in the case of mas- 
sive particles, stationary four- velocities. These four-velocity forms are then 
universally applicable on the given geometry (principle of equivalence), and 
can be used to properly incorporate geometrodynamics (kinematics) with- 
out additional entanglements being introduced to any (internal) quantum 
dynamics. The trajectories of classical particles were plotted and found to 
satisfy the constraints expected from causality on the Penrose diagram. 

In addition, simple quantum behaviors on the black hole background 
were seen to follow intuitively consistent behaviors. First, the form of the 
equation describing the dynamics of spherically symmetric solutions to the 
massless Klein-Gordon equation were developed. Numerical solutions of 
gravitating quanta on the black hole background were found to have den- 
sities most prevalent in the regions near the singularity and horizon. The 
numerically explored Klein-Gordon wavefunction was continuous through- 
out the global space-time, but its derivative was found to behave in a 
singular manner near the horizon and the ingoing horizon. The behavior 
of the energy-momentum density of the scalar field was then compared 
with that necessary to generate the dynamic black hole. It was found that 
the non-linear dependency of the Klein-Gordon energy-momentum tensor 



on field derivatives make it problematic for this scalar field to be used to 
consistently mimic the radial dependency of the Einstein tensor describing 
the dynamic black hole. 

Next, generic complex scalar fields whose phase relationships are lin- 
early related to energy-momentum densities were developed. This means 
that the energy-momentum relationships that generate coherent phase in- 
terference properties (as experimentally observed) directly relate to the 
energy-momentum tensor calculated for that field. The generic form can 
be generalized to spinor fields in a straightforward manner. The scalar 
form was shown to satisfy the expected physical dynamics of quantum ob- 
jects (e.g. probability flux conservation, phase coherence, etc.) on the 
black hole background. Numerical solutions for the quantum state vector 
describing the scalar fields depend only upon the boundary conditions of 
any given quantum state, and are straightforward to calculate. The regions 
of the global space time available for outgoing quanta are delineated by the 
black hole horizon, while regions available for ingoing quanta are delineated 
by the ingoing horizon of the black hole. The energy-momentum tensors 
of any of these quantum states linearly contribute to the overall system 
kinematics and dynamics. 

The construction of mutual/self gravitating macroscopic collections of 
the scalar fields was next explored. A straightforward development of co- 
gravitating solutions was expected to require the following properties: (i) 
cluster decomposability of non-interacting radiations/quanta, (ii) proper 
time/affine dynamics for substantive kinematic flows in order to incor- 
porate the principles of relativity/equivalence, and (iii) linear forms for 
additivity of the energy- momentum contributions. The solutions found in- 
volved the collective contributions of the energy-momentum densities of 
the mutually gravitating quanta combined with a real, spatially coherent 
core gravitating field resulting in a self-consistent gravitating system. The 
real core field for the excreting black hole examined consists of outgoing 
massless quanta which have a direct interaction term with the metric form 
in the core Lagrangian. Because of the direct metric interaction, the core 
fields do not have conserved probability densities, but have a form which is 
fixed by that of the Einstein tensor. The key points of flexibility that allows 
generic solutions are that the Euler-Lagrange equations for the core gravi- 
tating field only require that the field be massless, while the co-gravitating 
decomposable fields contribute to the overall energy- momentum densities 



in an additive way directly related to the quantum phase coherence ener- 
getics. 

Since the approach here taken is generic and not peculiar to the black 
hole examined, it should be capable of describing several physically relevant 
gravitational scenarios. In particular, this generic form should be useful 
for developing mutually gravitating solutions for static backgrounds in the 
exterior region (e.g. planets, stars, etc.) using ingoing and outgoing quan- 
tum stationary states, along with any kinematically stationary solutions. 
Fields (such as the generic gravitating quanta or gravitating clusters) with 
vanishing extremal Lagrangian can always be algebraically added to static 
geometries to help generate vacuum solutions to Einstein's equation. 

The author is involved with several related on-going explorations of 
the geometrodynamics of gravitating quanta. One might gain insight into 
the microscopic contributions of the individual quantum clusters by more 
closely examining the single quantum limiting forms of the macroscopic 
system (whether the system has a singularity or not). In addition, the 
exploration of the co-gravitating initiation of accretion to form a singularity 
might answer fundamental questions on how black holes start the process 
of collapse. Finally, the construction of mutually gravitating cosmological 
systems can be likewise accomplished in a straightforward manner, as will 
be presented shortly in a follow up paper. 
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